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Abstract 

The distribution of Ge-O-Ge and Si-O-Si bond angles a in amorphous germania 
and silica is re-determined on the basis of diffraction experiments. The bond angle a 
' joining adjacent tetrahedra is the central parameter of any continuous random network 

description (CRN) of these glasses. New high energy photon diffraction experiments 
on amorphous germania (at photon energies of 97 and 149 keV) are presented, cover- 
ing the momentum transfer 0.6-33.5A" 1 . In photon diffraction experiments on Ge02 
the contribution of the OO pairs is very small. To obtain a similar information for 
amorphous Si02, high energy photon diffraction experiments jlj have been combined 
with neutron diffraction data [S, B on amorphous silica in order to eliminate the OO- 



o 
m 
o 

On 



CD 

o 



- 1—1 

X 



partial structure factor. With this technique it is shown that the Si-O-Si angle dis- 
tribution is fairly narrow (a = 7.5°) and in fact comparable in width to the Ge-O-Ge 
angle distribution (a = 8.3°), a result which differs from current opinion. The narrower 
distribution found in this study are in much better agreement to the determinations 
based on 29 Si-MAS-NMR. Among the various models relating the chemical shift to the 
bond angle, best agreement is found with those models based on the secant model. 
Sharp components in the bond angle distribution can be excluded within the reached 
real space resolution of 0.09A. 



1 Introduction 

Vitreous silica and germania are two prototype simple glasses. These two glasses consist 
of corner-sharing tetrahedra linked to a continuous three-dimensional network. Numer- 
ous diffraction experiments have been published, especially concerning amorphous silica, of 
which the eldest date back to the very beginning of x-ray diffraction Q. But even before, 
some ideas about the structure of glass already existed 0] . The traditional theories of glass 
structure are the crystallite theory and the continuous random network (CRN) theory at- 
tributed to Lebedev |(| and Zachariasen [[tJ respectively. Modern crystallite theory predicts 
fluctuations of the intermediate range order where some regions approach the atomic ar- 
rangement of the crystal. Fluctuations of the intermediate range order are also allowed in 
the CRN model, so the difference between both models is mainly a gradual one. As bond 
angles in crystals show only small variations around discrete values, modern crystallite the- 
ory should lead to sharp bond angle distributions superimposed on the broader distribution 
prevailing in the more disordered interconnecting regions. A more detailed discussion of 
these topics is given in the review articles of Wright || and Porai-Koshits Q . 

The central parameter of a CRN model of germania and silica is the distribution of A- 
O-A bond angles (A=Si,Ge) joining adjacent tetrahedra. The Si-O-Si angle distribution in 
amorphous Si02 is often believed to be much broader than the corresponding distribution 
in germania. This can be traced back to the classical work of Mozzi and Warren . Later 



some discussion came up about the average Si-O-Si angle |Tl|, |T^|, but not about the width 
of the distribution. As the mean Ge-O-Ge angle of 133° is close to the one for planar 3- 
membered rings, Barrio et al. |]l3f recently argued, that the narrower distribution in Ge02 
is an indication for the presence of such small, stiff ring-systems in a much larger extent 
than in Si02- This hypothesis was inspired by the interpretation of IR-data of these glasses. 
Robertson et al. [jwj computer relaxed the Bell and Dean CRN model, but were quite 
suspicious on their result, as the Si-O-Si angle distribution found was much narrower than 
the Mozzi and Warren distribution. 

The Si-O-Si angle can also be determined with 29 Si-MAS-NMR spectroscopy. The 
method is based on the assumption that a functional relation exists between the chemical 
shift 8 nmr and the Si-O-Si bond angle. This relation has to be established experimentally in 
measuring the chemical shifts of a sufficiently large number of crystalline Si02 polymorphs 
with known Si-O-Si bond angles. The first Si-O-Si angle determination by this technique 
was described by Dupree and Pettifer jlq |. This first distribution was claimed to be very 
broad (FWHM 40° !) and essentially flat between 140-155°. This result was later on revised, 
being attributed to boron, sodium and aluminum contaminants present in the investigated 
glass sample jT3|. In this paper also the influence of various assumptions on the relation 
between the chemical shift and the bond angle is discussed. In general, the bond angle 
distributions determined with NMR are found to be substantially narrower 0, H| than the 
Mozzi and Warren distribution, although in jl8| a good agreement with Mozzi and Warrens 
result is claimed (This is certainly to be understood in comparison to flq] ). MD simula- 
tion sometimes agree with the NMR results |Tgf| , while other recent simulations |Q find 
good agreement with Mozzi and Warrens result and a third group reported even a broader 
distribution [pl| . 

Recently, diffraction experiments on amorphous systems (glasses: 0, molecular liquids: 
p^|) with high energy photons as delivered from modern synchrotron radiation sources 
became feasible. The term high energy photons will be used for photons in the energy 
range of ~ 80 — 200 keV, i.e. of an energy about one order of magnitude higher than 
that of x-rays traditionally used in angle dispersive measurements. As the photo-electrical 
absorption decreases with energy approximately like E 3 , the absorption is reduced by three 
orders of magnitude. Scattering therefore becomes the dominant process in a large number 
of materials and the conditions of the experiment become in general quite similar to that 
of a typical neutron diffraction experiment. The wavelength of the photons is quite short 
and consequently no principle limit exist on the maximum momentum transfer attainable: 
The momentum transfer for back scattering with a photon energy of 100 keV is Q ma x = 
4%E/(hc) = AttE ■ 0.080658keV _1 A 1 = 101. 4A~\ Considering only a more commonly 
used momentum transfer range the scattering takes place through small or intermediate 
scattering angles. This is advantageous for the use of small windows or flat area detectors, 
but the first two characteristics - low absorption and high momentum transfer - will be 
more important for the present study. For a description of the use of high energy photons 
in general - i. e. not restricted to amorphous substances - see f23f . 

This article will focus on the Si-O-Si and the Ge-O-Ge bond angle distributions and its 
comparison. The data on amorphous germania are new and the experiment is described in 
the following section. The experiment on vitreous silica at the DORIS-III synchrotron has 
been described previously 0. The neutron diffraction data on amorphous silica have been 
taken from Q and ||. 
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2 Experiment on ger mania 



The vitreous germania has been prepared from 99.999% pure GeC>2 (Strem- Chemicals, New- 
buryport, MA, USA) by melting the crystalline material. The melt was kept for 8 h at 1473 
K at atmospheric pressure (melting point 1383 K). The furnace was then switched off and 
the sample cooled to 900 K in about 10 min. The resulting sample was clouded as a result 
of the inclusion of small gas bubbles but no traces of crystallinity were left. 

The diffraction experiment on vitreous germania has been carried out at the high energy 
beam-line ID15A of the European Synchrotron Radiation Facility ESRF. The Synchrotron 
was operated in 16 bunch mode with an average current of about 70 mA and beam life- 
times of about 16 h. The beam-line is designed to operate alternatively with a permanent 
magnet multipole wiggler or with a super conducting wavelength shifter ]24| . The former 
only is commissioned at the moment and thus has been used for the experiment. Its gap 
was maximal closed to 20.3 mm leading to a critical energy of 43 keV. To remove the heat 
load, lower energies are cut by a cooled 4 mm aluminum filter. The beam size can be chosen 
by a set of water-cooled secondary slits in front of the monochromator position in the optics 
hutch II 0. 

A perfect (220) Si crystal in Laue geometry was used for monochromatization and 
mounted on a bending device. To reduce background in the experimental hutch, the white 
beam was stopped in the optics hutch II. 

A dedicated two axes diffractometer at the high energy beam-line does not exist. Instead 
the sample was mounted on the analyzer stage of the triple axes diffractometer. Here, 
the scattering angle is chosen by a translation of the detector perpendicular to the white 
beam axis. This geometry allows to select a maximum scattering angle of 29 = 26°. The 
solid state detector was equipped with a 30 mm thick germanium crystal and has a typical 
energy resolution of about 1 keV. Discriminating the energy around the photon energy of 
the monochromatic beam reduces the background scattered from the white beam. The 
intensity in the first sharp diffraction peak of the sample was around 5 • 10 4 counts per 
second. Further increase of the intensity would raise considerable dead-time problems. In 
order to cover a large momentum transfer range two different photon energies of 97.58 keV 
and 149.80 keV were chosen. The monochromator was bent for the higher energy setup in 
order to compensate for the intensity drop. 

The energy of the primary beam was determined by observing the Bragg peak position of 
a second perfect (220) Si crystal mounted on the monochromator position of the triple-axes 
diffractometer for both, dispersive and non-dispersive geometry. 

Due to the excellent beam stability of the ESRF synchrotron, normalizing the data 
on the beam current was found to be more convenient than using an Nal monitor, which 
was installed for control purposes. The latter observes the air scattering produced by the 
incoming beam and was too sensitive to the hutch background as its energy resolution is too 
poor. The problems associated with beam fluctuations described in jjj are mostly reduced 
in the present study. In order to minimize generally the effect of drifts the data acquisition 
was organized as a sequence of several short scans. 

The data were subsequently corrected for background, absorption, multiple scattering, 
polarization of the incoming beam, the variation of the solid angle seen by the detector as a 
result of the changing distance between the sample and the detector (tangential movement, 
see above) and the increased detection efficiency for inelastically scattered photons, [ir 
is 0.25 at 149 keV (with n being the total absorption coefficient and r the radius of the 
sample). Air scattering is the dominant background contribution, especially at low Q's. 
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The correction procedure was described in detail in The resulting fully corrected and 
normalized data are presented in Fig. 1. 

3 Data analysis 

The usual definition of the structure functions in the neutron and photon case are: 

da\ n ( da\ x 

c,,,» - ^ dQ )d,st y bjbk \dn) dlst \- fjfk (1) 
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where t is the distinct (interference) differential cross section, S and bi are the struc- 

ture function and the scattering lengths in the neutron case, i and fi are the structure 
function and the form-factors in the photon case, the sum is extending over the unit com- 
position. S(Q) and i{Q) are weighted sums of the partial structure factors Sj k - The Sjk 
are, of course, the same in neutron and photon diffraction. The weighting coefficients in the 
photon case are constant in Q only if one assumes that the form factors of the atoms present 
differ by a scaling factor only - this is the K-approximation introduced by Warren, Krutter 
and Morningstar [ p5| . A combination of neutron and photon data given by: 

(da \ n ( da \ x 

— J ~ V & fc' ( 3F7 ) = ifj'fk'bjbk-bfb k >fjfk)s jk (2) 

ai ' / dist \aiij dist fe( 

eliminates one particular partial structure factor Sfk 1 from the diffraction pattern without 
any approximation. This procedure is analogous to the isotope substitution technique in 
neutron diffraction, but has the advantage that the same sample can be used. Prerequisites 
for a successful separation are a good knowledge of the form-factors and of the neutron 
scattering lengths. The neutron scattering lengths used in the analysis are taken from p(| , 
the form- factors from |27], |2§|] . 

The three differences which are possible for SiC>2 according to Eq. (0) arc shown in Fig. 
2. The contribution of the first two peaks in real space - the SiO- and 00- first shell - 
calculated according to || is indicated. It is apparent from the figure that the scattering 
pattern is dominated by these peaks above 15 A -1 . It is also apparent that a systematic 
error is present between 16 A -1 and 20 A -1 , probably caused by the splicing of two scans 
in that region Q. This can also be seen in Fig. 3 where the contribution of the first two 
peaks has been subtracted. 

A useful real space correlation function which can be used for the analysis of amorphous 
materials is (in the neutron case): 

T(r) = -j^ J^* QS(Q)M(Q) sin(Qr)dQ + r = ^ -g(j^t' jk (r), (3) 

where M(Q) is an optional modification function. T(r) is related to g(r) - the correlation 
function which is more commonly used in the discussion of the structure of liquids - simply 
by T(r) — rg(r). The experimental component correlation functions t' are related to the 
real component correlation functions t via a convolution with a peak shape function [[l0| , 
p9| . The peak shape function is determined by the modification function used and - in 
photon diffraction - the non-constant weighting factors of the partial structure factors. The 
modification function used here is the rectangular function (all data points are equally 
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weighted for the Fourier transform): 

M(Q) = 



Q ^ Qmin 
Qmin ^ Q ^ Qr, 

Q> Q max 



(4) 



This is to the authors opinion a better solution than weighting parts of the measured data 
with a modification function as e. g. the one proposed by Lorch ]3Cj as long as the implica- 
tions on T{R) are noticed. 

For the differences given by Eq. (||) the Fourier transform has been calculated according 

to: 



T d iff(r) = 



1 



Wptfj'fubjnbu, - b j >b k >f j »f k ,,)N 



Qdj.v (Q)M(Q) sin(Qr)dQ + r (5) 



The division by fj'fk'bj"bk" — bjibk* fj» fk" removes the dependence of the peak shape func- 
tion from the non-constant weighting factors for one of the partial correlation functions (for 
fk"). N is defined by: 



/ J v(0)/ fc ,(0)6 i 6 fc -6 J v6v/i(0)/ k (0) 



(6) 



Whenever comparisons between a model and the experiments are presented, a double 
Fourier-transform has been carried out: The model partial correlation functions t(r) are 
transformed into Q-space: 



poo 

Qs{Q) = 4tt P / t(r) sin(Qr)dr 
Jo 



(7) 



multiplied by the appropriate weighting factors, summed and back transformed to r-space 
with the same Q-interval than the experimental data. In short, the model correlation func- 
tions are convolved with the appropriate peak shape function before comparison. The agree- 
ment with the experiment is judged according to || by calculating: 



frZ:: Tl p (r)dr J rmt 



(8) 



Models have been developed using a MC algorithm. The A-0 bond length (with A=Si,Ge) 
is assumed to have a Gaussian form: 



K ■ exp 



{ta-o - rp) 2 
2ol 



(9) 



with K being a scaling factor, ta-o the bond length, r Q and cr r having its usual significance. 
The bond angles A-O-A and O-A-0 have been modeled by distributions of the form: 



V(a) — K ■ exp 



(a - a ) 2 
2a 2 



■ sin(a), 



(10) 



with K again a scaling factor, a the A-O-A bond angle, ao and a are the free parameters of 
the distribution. The O-A-0 bond angle will be called (3 in the following, and a has to be 
replaced by (3 then in the above equation. The dihedral angle distribution has been given 
the form: 

[S - S a ■ (i - 2)p 



V(8) = K ■ exp 



2*1 



(11) 
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where <5 is the dihedral angle A-O-A-0, So has been set to 60°, K is a scaling factor and as is 
a variable parameter. This leads to distributions having the necessary threefold symmetry 
favoring eclipsed conformation (In the limit as —> oo an uniform dihedral angle distribution 
is reached). Assuming independence of the above distributions the short to intermediate 
range order is specified: The A-O, O-O, A-A first and second shell contributions. These 
distributions are calculated and compared with the experiment via the above described 
formalism. 



4 Simulation 

To calculate the distance distributions from the distributions (^) -(|TT|) chains O-A-O-A-O-A 
of A and O atoms were generated. Beginning the chain from one side - say with the O- 
atom a (normally distributed) random number determines the position of the next A- atom. 
The distribution (^|) is, hence, fulfilled in this step. The position of the next O- atom is 
determined by two random numbers (representing distribution ([)]) and (^)), the following 
A- atom by three random numbers and so on. As only four- body correlations and no cross 
correlations are considered a maximum of three random numbers is needed for positioning 
an atom in the chain. Each chain yields five times a A-0 first shell distance, two times 
a 0-0 and A-A first shell, three times a A-0 second shell and one times a 0-0 and A-A 
second shell distance. 10 5 chains were generated for each model. 

This MC approach has the advantage over the analytical approach followed in Q| to 
allow a much more flexible choice of distributions, while still being computationally very 
cheap (few seconds of CPU time) . It has the disadvantage that the analytical form of the 
distributions has to be imposed in fore-hand. It is easy to overcome this problem, following 
the ideas of reverse Monte- Carlo (RMC) simulations |3l[]. First one generates chains with 
quite wide distributions of next neighbor distances, bond angles and dihedral angles. The 
generated chains are accepted or rejected using the Metropolis criterion, if these chains 
improve or deteriorate the agreement of the total model (the ensemble of chains accepted so 
far) with the experiment. This technique is by orders of magnitude more time consuming as 
it needs the double Fourier transforms described in the preceding chapter to be executed for 
each chain generated. The next more complex simulation would be a full RMC simulation 
of the system: This introduces an additional information into the simulation process , i. 
e. the macroscopic density, but on the other hand - a disadvantage in a system which 
is in reality of course aperiodic - periodic boundary conditions. A RMC simulation for 
Si02 has already been carried out (including information from both neutron and photon 



diffraction experiments) 32 but it leads to structures which are apparently unphysical on 
various counts. One reason is perhaps, that a "good" structure can not be reached within 
the computer time available. 



5 Results and Discussion 

T x {r) for amorphous germania are shown in figure 4a. Looking closer at the GeO next 
neighbor peak in Fig. 4b, one recognizes that the peak width decreases and the peak 
height increases steadily with increasing Q ma x ■ The peak form is largely determined by the 
truncation effect and the Fig. 4b gives, hence, an idea of the real space resolution of various 
techniques: Q ma x — 14A corresponds to a conventional x-ray diffraction experiment, 
Qmax = 23A to the D4b neutron diffractometer of the ILL (Grenoble) and finally Q m ax = 
33 A can be easily reached with high energy photons. This behavior also proves, that the 
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region 23A < Q < 33A still contains valuable information about this peak. The GeGe 
next neighbor distance in Fig. 4c, in contrary, does not change its form for Qmax > 23A _1 . 
The GeGe peak has consequently reached its natural form and extending the integration to 
higher Q values can only obscure the peak shape by inclusion of the high Q noise. From 
the position and the width of the GeO and the GeGe next neighbor peaks one can deduce 
the bond angle distribution. o.q = 133.0° and a = 8.3°. This is in good agreement with 
an earlier x-ray structure determination fl33| using Mo-K" x-rays (17.4 keV), ending with 
the more qualitative statement: "[...] the average Ge-O-Ge angle is 133°. The distribution 
of Ge-O-Ge angles is fairly narrow but cannot be determined quantitatively with data of the 
present resolution." \ defined by Eq. (Q) is 2.0% in the range 2.5A < r < 4.0A and 2.8% 
in the range from < r < 4A. \ 2 does not change significantly, when the o~s parameter is 
changed from infinity (randomly distributed dihedral angles) to 38°. Decreasing as further 
deteriorates the fit. In the limit of pure staggered or pure eclipsed conformation (as — » 0) 
a x of 6.5% and 5.9% is reached. The final model parameters are collected in Tab. 0. The 
model is compared to the experimental Qi(Q) in Fig. 5. 

There is no indication for a sharp component in the bond angle distribution, which should 
be expected, if a significant volume fraction of the glass would reach crystalline order. 

The question arises, whether the connection of the tetrahedral units is really that different 
in germania and silica. To answer this question a reanalysis of previous high energy photon 
experiments jjj was started, combining them with neutron diffraction data (2[ |^| in the way 
described above. The most interesting of the first order differences dyw presented above is 
the one where the 00 contribution is removed. The Fourier-Bessel transform is shown in 
Fig. 6. Qmax has been chosen to be 16A" 1 for the reasons outlined in the previous section. 
The data are normalized according to Eq. ^ such that the SiSi partial structure factor has 
a constant weighting factor. Assuming again Si-O-Si bond angle distributions as in Eq. [l0|, 
broad distributions with a > 10° are clearly unacceptable. From that it follows that the 
assumption of isotropic dihedral angles is not justified. The geometric relations are such, 
that isotropic dihedral angles and narrow bond angle distributions lead to a shoulder or even 
a distinct peak on the low r side of the SL..0 second neighbor peak. Some evidence was 
already collected by earlier authors on the non-randomness of the dihedral angles, but the 
literature is controversial in this point. So, c. g. Galccncr |34| (from sterical arguments) and 
Evans et al. (from a comparison of various models with neutron diffraction data) suspected 
the dihedral angles to be non-randomly distributed. Gaskell et al [ |36] found a non uniform 
dihedral angle distribution in his relaxed version of the Bell and Dean pl| model, while 
Gladden (3^] on the other hand found uniform dihedral angles in her models. The best 
agreement with the experiment is reached with ao = 148.3°, a — 7.5° and <7,5=27°. The 
final parameters of the model are collected in Tab. 1 

The interpretation of the NMR line-shape of silica glass is based on the knowledge of a 
relation between the chemical shift and the Si-O-Si angle. Four such relations are discussed 
in reference 0: The secant Q, the linear |3^] ,the point charge Q and the s-hybridization 
]4l| | model. The secant and the s-hybridization model produce almost identical Si-O-Si angle 
distributions and are not considered separately. The resulting bond angle distributions are 
compared in Fig. 7 with each other, and with the distribution obtained in this work. \ 2 as 
a function of r max is shown in Fig. 6. While the linear model gives a bad agreement with 
the combined neutron and photon data, the secant (with the data of Pettifer et al.) and 
point charge model give a comparable fit quality , but the lowest \ °f au NMR distributions 
has the bond angle distribution given by Gladden et al. These authors also used the secant 
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model to derive their bond angle distribution. It is interesting to note, that the tail to 
large angles present in many published Si-O-Si distributions is not necessary to explain the 
diffraction data. The x values reached by various distributions are reported in Tab. 2. All 
experimenters, of course, used different samples - with different concentrations of impurities 
and defects, different quenching conditions etc. - and this may also play a significant role 
at this level. It would be interesting to go an alternative route: The determination of the 
functional relation between chemical shift and bond angle by investigating the same glass 
sample with neutron and photon scattering and 29 Si- MAS- NMR. 

Several authors devoted special attention to the first sharp diffraction peak (FSDP) (see 
the discussion given by Elliot It is noted, that the FSDP vanishes almost completely 

when the SiO partial structure function is eliminated (Fig. 2). Probably the 00- and the 
SiSi- partial cancel as its weighting factors have different sign. Furthermore it is noted that 
the FSDP can not even be recognized any more if the contribution of the short range order 
is subtracted (Fig. 3). It is thus concluded that the FSDP by itself has very little structural 
significance as it is more or less arbitrarily composed of the weighted partial distribution 
functions and contributions of short and medium range order. This will be even worse in 
glasses of more complicated composition. 

The experimental QS(Q) and Qi(Q) respectively are compared to the model in Fig. 
8. The specified elements of the short and intermediate range order can describe both the 
neutron and photon diffraction pattern down to about 2.5 A~ x . The contribution of the 
different first and second neighbor distances to T x (r) is shown in Fig. 9 and can be compared 
to Mozzi and Warrens Fig. 4 fl(i|| . 

Finally, one should discuss what possibilities exist for a further disentangling of the 
partial structure factors. There is very little that can be done for silica. There are no 
suitable isotopes for oxygen, the scattering length difference between 28 Si and 29 Si is small 
(0.6fm) and 29 Si is quite expensive. Anomalous diffraction and EXAFS spectroscopy are 
hindered by the low energy of the Si-K edge (1.8 keV). This is, why a combination of photon 
and neutron data is especially valuable in this case. In contrast, there are suitable Ge- 
isotopes for an isotope substitution experiment ( 70 Ge and 73 Ge) [^6|. So, Ge02 could be a 
promising candidate to get a very detailed structural picture of a simple tetrahedral glass. 

6 Conclusion 

Combining neutron and photon diffraction can yield additional information on such well 
studied systems as the simple glasses. The ideal counterpart to modern neutron instru- 
mentation are high energy photons in the energy range from ~ 80-200 keV, which can be 
delivered from modern synchrotron radiation sources. The main advantages which come 
into play in the present study is the large accessible Q-range and the low absorption. 

Due to the large accessible Q-range the GeGe next neighbor peak can be observed without 
any truncation effect. This leads to an accurate Ge-O-Ge bond angle distribution, the central 
parameter of a CRN description of the glass. There is no indication of a sharp component 
in the bond angle distribution, which could be expected, if a significant volume fraction of 
the glass would reach crystalline order. 

Combining neutron and photon results can lead to an isolation of the SiSi first neighbor 
peak in Si02, too. From this, it can be deduced that the very broad Si-O-Si distribution 
calculated by Mozzi and Warren is not acceptable. This broad distribution is a consequence 
of the assumption of randomly distributed dihedral angles. Likewise the photon data used 
here would lead to a broad Si-O-Si bond angle distribution - when interpreted alone without 
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combination with neutron data and using this wrong assumption Only in combination 
with neutron data, it becomes clear that preference should be given to narrower distribu- 
tions and the width of the Si-O-Si distribution is indeed within the experimental accuracy 
found to be identical to the width of the Ge-O-Ge angle distribution. From the bond angle 
distributions determined by 29 Si-MAS-NMR the one given by Gladden et al. gave the best 
fit to the combined neutron and photon data. The Si-O-Si distributions determined form 
diffraction and from NMR measurements are in much better agreement, when neutron and 
photon diffraction are combined. Among the models relating the chemical shift to the bond 
angle the secant model agrees best with the diffraction data. These narrower distributions 
demand in turn a non-uniform distribution of the dihedral angles, that is to say a preference 
of the staggered conformation. The bond length, the Si-O-Si angle and the dihedral angle 
are likely to be correlated among each other and in neighboring units, but no evidence can 
be taken from diffraction experiments alone on these correlations. 



7 Acknowledgments 

Financial support of the Deutsche Forschungsgemeinschaft DFG under Grant No. Ne584/1-1 
is gratefully appreciated. 



References 

[1] H. F. Poulsen, J. Neuefeind, H.-B. Neumann, J. R. Schneider and M. D. Zeidler, J. 
Non-Cryst. Sol. 188 (1995) 63 

[2] P. A. V. Johnson, A. C. Wright and R. N. Sinclair, J. Non-Cryst. Sol. 58 (1983) 109 

[3] D. I. Grimley, A. C. Wright and R. N. Sinclair, J. Non-Cryst. Sol. 119 (1990) 49 

[4] P. Debye and P. Scherrer, Nachr. d. Konigl. Ges. d. Wiss. z. Gottingen, Math.-phys. 
Klasse (1916), 16; S. Kyropoulos, Z. anorg. allg. Chcm. 99 (1917) 197 

[5] J. Tyndall,"Heat- A Mode of Motion.", London(1863); M. L. Frankenheim, "Die Lehre 
von der Cohasion", Breslau (1835), 389 

[6] A. A. Lebedev, Trudy Gos. Opt. Inst. 2 (1921) 1 

[7] W. H. Zachariasen, J. Am. Chem. Soc. 54 (1932) 3841 

[8] A. C. Wright, PAC RIM Meeting, Honolulu, HI, USA, Nov. 1993, J. Non-Cryst. Sol. 
179 (1994) 84 

[9] E. A. Porai-Koshits, J. Non-Cryst. Sol. 123 (1990) 1 

[10] R. L. Mozzi and B. E. Warren, J. Appl. Crystallogr. 2 (1969) 164 

[11] J. R. G. da Silva, D. G. Pinatti, C. E. Anderson and M. L. Rudee, Philos. Mag. 31 
(1975) 713 

[12] P. G. Coombs, J. F. Natale, P. J. Hood, D. K. McElfresh, R. S. Wortman and J. F. 
Shackelford, Phil. Mag. B51 (1985) L39 

[13] R. A. Barrio, F. L. Galeener, E. Martinez and R. J. Elliott, Phys. Rev. B 48 (1993) 
15672 



9 



[14] J. L. Robertson and S. C. Moss, J. Non-Cryst. Sol. 106 (1988) 330 

[15] R. J. Bell and P. Dean, Phil. Mag. 25 (1972) 1381 

[16] E. Dupree and R. F. Pettifer, Nature 308 (1984), 523 

[17] R. F. Pettifer, E. Dupree, I. Farnan and U. Sternberg, J. Non-Cryst. Sol. 106 (1988) 
408 

[18] L. F. Gladden, T. A. Carpenter and S. R. Elliot, Phil. Mag. 53 (1986) L81 

[19] P. Vashishta, R. K. Kalia, J. P. Rino and I. Ebbsjo, Phys. Rev. B 41 (1990) 12197 

[20] B. Vessal, M. Armini and C. R. A. Catlow, J. Non-Cryst. Sol. 159 (1993) 184 

[21] B. P. Feuston and S. H. Garofalini, J. Chem. Phys. 89 (1988) 5818 

[22] J. Neuefeind, M. D. Zeidler and H. F. Poulsen, Mol. Phys. 87 (1996) 189 

[23] J. R. Schneider, Synchr. Rad. News 2(8) (1995) 26 

[24] ESRF beam-line handbook Dec. 94, ed. R. Mason, p. 87; World Wide Web: 
http: / /www. esrf.fr 

[25] B. E. Warren, H. Krutter and O. Morningstar, J. Am. Ceram. Soc. 19 (1936) 202 

[26] L. Koester, H. Rauch and E. Seymann, At. Data Nucl. Data Tabl. 49 (1991) 65 

[27] J. H. Hubbell, W. J. Veigele, E. A. Briggs, R. T. Brown, D. T. Cromer and R. J. 
Howerton, J. Phys. Chem. Ref. Data 4 (1974) 471 

[28] J. H. Hubbell and I. Overb0, J. Phys. Chem. Ref. Data 8 (1979) 69 

[29] J. Waser and V. Schomaker, Rev. Mod. Phys. 25 (1953) 671 

[30] E. A. Lorch, Phys. Rev. C 2 (1969) 229 

[31] R. L. Mc Greevy and L. Puztai, Molec. Simul. 1 (1988) 359 

[32] D. A. Keen and R. L. McGreevy, Nature 344 (1990) 423 

[33] A. J. Leadbetter and A. C. Wright, J. Non-Cryst. Sol. 7 (1972) 37 

[34] F. L. Galeener Phil. Mag. Letters B51 (1985) LI 

[35] K. M. Evans, P. H. Gaskell and C. M. M. Nex in 'The Structure of Non-Crystalline 
Materials, eds. P. H. Gaskell, J. M. Parker and E. A. Davis (1983), p. 426 

[36] P. H. Gaskell and I. D Tarrant, Phil. Mag. B 42 (1980) 265 

[37] L. F. Gladden, J. Non-Cryst. Sol. 119 (1990) 318 

[38] J. V. Smith and C. S. Blackwell, Nature 303 (1983) 223 

[39] G. Engelhardt and R. Radeglia, Chem. Phys. Lett. 108 (1984) 271 

[40] J. M. Thomas, J. Klinowski, S. Ramdas, B. K. Hunter and D. T. L. Tennakoon, Chem. 
Phys. Lett. 102 (1983) 158 

[41] U. Steinberg, Mol. Phys. 63 (1988) 249 

10 



[42] S. R. Elliot, Nature 354 (1991) 445 

[43] R. N. Sinclair, D. A. G. Johnson, J. C. Dorc, J. H. Clarke and A. C. Wright, Nucl. Inst. 
Meth. 117 (1974) 79 



11 



8 Figures 



Figure 1: The fully corrected and normalized data (high energy photons, Ge02) 

Full line: experimental intensity; dashed line: isotropic part of the scattering intensity 

(Self scattering + Compton scattering) 

The insert is a zoom at the high-Q region. 

Figure 2: The first order differences dfk> in amorphous silica 

From top to bottom: QdsiSi(Q), Qdoo{Q), QdsioiQ)- Solidline: experimental Qdj>k> (Q), 
dashed line contribution of the first two peaks (SiO-,00-) in real space. 

Figure 3: The first order differences in amorphous silica, contributions of the first two peaks 
in real space removed 

From top to bottom: Qd Sl Si{Q), QdooiQ), Qdsio{Q)- 

Figure 4: T x (r) for vitreous germania 

a) T x {r) for vitreous germania, Q m ax — 23A 

b) and c) Closer look at the GeO- and the GeGe-peak respectively: solid line Q m ax = 
33A , dashed line Q ma x = 23A , dots Q ma x = 14A 

Figure 5: Comparison of the model distribution for Ge02 with experiment 
Solid line: experimental Qi(Q) , dashed line: model Qi(Q) 

Figure 6: Too( r ) and X 2 horn the various distributions 

Upper window:Too i r ) > lower window x 2 {\ 2 is not normalized to j T 2 xp (r) here). 

Full line: experiment, dashed triple dotted: Mozzi and Warrens distribution p0|] , 1. 
dotted line: Best fit of Eq. (10), 2. dashed dotted line: NMR model of ref . @, 3-5: NMR 
models discussed in ref. |l7j]; 3. linear model (not shown) 4. dashed double dotted line: 
NMR data of ref. jL7), point charge model; 5. dashed line NMR data of ref. fll7fl , secant 
model; 6. triple dotted line: best fit of Eq. (10) to the NMR distribution of ref. Q 

Figure 7: Various models for the Si-O-Si bond angle 

1. dotted line: Best fit of Eq. (10) to the combined photon and neutron data, 2. dashed 
dotted line: NMR result from ref. |l8), based on the secant model; 3-5: NMR models 
discussed in ref. p!7| ; 3. dashed triple dotted line: linear model, 4. dashed double dotted 
line: point charge model, 5. dashed line: secant model. 6. triple dotted line: best fit of Eq. 
(10) to the NMR result from ref. Q 

Figure 8: Comparison of the model for Si02 with the experiment 

Solid line: experimental QS(Q) and Qi(Q) respectively; dashed line: model QS(Q) and 
Qi(Q) respectively, a) photon diffraction jlj, b) neutron diffraction ||, |j| 

The neutron data were available as spline fit to the experimental data only. 

Figure 9: First and second shell contributions to T x {r) 

Solid line: experimental T x (r) and Fourier transform of the residue with Qmax = 
10.6A" , dashed line: SiO- and 00- first neighbor, dotted line: SiSi- first shell, dashed 
dotted: SiO- second shell, triple dotted line: OO second shell, double dotted line SiSi- 
second shell contribution 
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Table 1: Structural parameters for SiC>2 and GeC>2 





TAX 


a r 


ao cr a 


00 




crs 


Ge0 2 


1.73 


.0415 


133.0 8.3 


109.47 1 




38 


Si0 2 


1.605 3 


.0493 3 


148.3 7.5 


109.47 1 


4.2 3 


27 



Distances in A, angles in degree. 

1: Not refined, /?o is taken from the ideal tetrahedron 

2: Not refined, taken to be equal to SiC>2 

3: Not refined, taken from reference || 



Table 2: \ reached from various bond angle models 



This work 


0.84% 


Gladden et al.\l^ 


1.4% 


Pettifer et al. 17 


secant model 
point charge model 
linear model 


2.26% 
2.20% 
3.41% 
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